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Abstract 


A  class  of  data  windows  is  proposed  for  t»sc  in  the  estimation  of  the  power 
spectra  cf  stationary  stochastic  processes.  These  windows,  which  are  genera  11  na¬ 
tions  of  the  standard  Farzeu  filter,  m o  constructed  by  the  vise  of  appropriately 
normalized  B-splir.os.  It  is  demonstrated  how  this  window  e]n.;s  nay  be  computa¬ 
tionally  Impl  t*’:*ont  c»J  using  c  East  Fourier  Trans  Form  algorithm,  7  ho  efficiency  cf 
the*  resulting  ptocederc  is  general  Iv  a  significant  improvement  over  Lhe  state  of 
the  art  v<th  little  additional  computer  tine  required. 


1.  ISTRiWCrrON 

TM  5  psprr  fs  concerned  with  the  problem  of 
digitally  estimating  the  pewer  spectrum  of  a  i'i  de¬ 
fense  stationary,  ervodic,  stochastic  process 
from  a  sample  function  of  finite  length.  Not 
long  after  Cooley  and  Tuhey  [5]  introduced  the 
fast  Fourier  Irens  form  algorithm  (FFTl,  its  appli¬ 
cation  to  the  estimation  of  power  spectra  was  dis¬ 
cussed  by  Bingham,  Godfrey,  and  Tukey  [2j  and  by 
Welch  [18],  They  suggested  a  method  based  on  com¬ 
plex  demodulation  and  pointed  out  the  computational 
speed  advantages  which  can  he  obtained  by  using 
the  FFT  for  computing  Fourier  poriodograms .  Welch 
proposed  a  direct  method  which  has  become  the 
standard  numerical  appioach  to  the  problem  of 
power  spectral  estimation.  This  procedure  has 
three  parts:  subdividing  the  available  data,  using 
a  data  window  t*1  compute  a  smoothed  spectral  esti¬ 
mate  for  each  segment,  and  averaging  these  smoothed 
spectral  estimates.  The  present  paper  proposes 
a  class  of  da  *  a  windows  derived  from  the  B-spllne 
basis  functions  of  Schoenberg  (14).  Krrpi  r  lea  I  ly , 
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the  result  of  usiii^  Llase  Grtu  windows  is  c  sig¬ 
nificant  improvement  in  estimates  of  power  spectia 
without  a  noticeable  increase  in  the  commit  at  i  c»nal 
effort.  We  have  chosen  to  u*c  continuous  notation 
throughout;  conversion  to  discrete  formulae  where 
appropriate  is  straight  forward. 

2.  PRELIMINARIES 

We  first  outline  some  basic  notation  and  re¬ 
sults  which  are  assumed.  Here  we  follow  Stein  ard 
Weiss  (16).  I.1  denotes  the  Banach  space  of 

Lebesgue  measurable  functions  which  nrc  (Lobes guc) 
integrable  over  the  real  line  V\  .  For  ,  the 

Fourier  transform  of  f  is  the  function  defined 
by 

kc)  -  r„  f(t)c'2"ictdt 

for  all  F,  €  P  where  l  =■  /-I  ,  The  notation 
(...)A  denotes  the  Fourier  transform  of  (...). 

For  f^GL1  ,  the  convolution  h  -■  f*g  is  defined 
as  the  function  h  C  L1  given  by 
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MT>  -  f(--t)g(t)dt 

fur  .ill  ?  6  F.  If  f,gr  l.'  ,  then  (f*>;)  -  fg. 

Similarly,  (fg)^  -  f-*g  provided  Che  functions  in¬ 
volved  are  such  that  all  operations  make  sense. 

For  a  >0,  represents  a  dilation  operator, 

i.e,,  \f(t)  =  f(at).  Whenever  f,gy  Ll,  it 

is  clear  that 

<V*«.8>(n  =  l  M(f*g)(?)l 

and  0,f)*<?)  =  i^v  f(’)  . 

We  shall  consider  a  sample  function  x(t) 

T  T 

on  a  finite  interval,  -  j  <  t  <  j  unless 
otherwise  stated,  of  a  stochastic  process 
( X(  t )  :  -«  <  t  <  .  We  assume  this  underlying 

stochastic  process  to  be  wide-sense  stationary 
(Doob{8])  and  ergodic.  For  simplicity  we  take 
the  mean  value  to  be  zero,  bet  X(c)  have  the 
spectral  density  function  P(f).  The  terminology 
related  to  spectral  analysis  used  in  this  paper  is 
generally  that  of  Blackman  and  Tukey  *[  3 )  . 

it  is  apparent  that  high  resolution  and 
stability  (snail  bias  and  variance)  are  desired 
qualities  o f  puwe r  spec t  ra  l  ns t ima t os .  However , 
a  compromise  between  the  two  must  be  made.  In 
order  to  eh tain  a  stable  spectral  estimate  using 
the  classical  indirect  method,  an  appropriate 
function  of  lag  must  modify  the  sample  autcco- 
variancc  function.  Fundamentally  the  indirect 
method  estimates  P(f)  via 

e.  (f)  -  r  C(t)l.(T)e*WfT  dr 

t .  m  *  v  „  r 

whore  t.(T)  is  the  finitely  supported  lag  window 
and  C  ( ^  is  the  sample  autocovar  i  mice  function 
,  (T- It |)/2  . 

<:(T)  F-(T-It|)/2  X(t+  J1x(t'2  ),lt- 

Usually  the  support  of  I.(T)  is  much  less  than 
the  entire  interval  -T<  r<.  T  . 

The  at  .uments  for  using  lag  windows  (e.g., 
see  Jenkins  and  Watts  [10)3  ran  be  adapted  to 
sliow  that  a  p:occdure  which  utilizes  the  FFT  to 
estimate  power  spectra  directly  should  incorporate 
a  data  window.  That  is,  for  a  direct  procedure 
to  have  a  proper  balance  between  resolution  and 
stability,  the  data  itself  must  be  modified  by  a 
■ultablr  function  of  time.  Thus,  basically  the 


at  discrete  frequencies  where  W(t)  is  the 
finitely  supported  data  window. 

The  end  result  of  either  of  these  modifica¬ 
tions  is  that  we  in  fact  estimate  smoothed  values 
of  P(f).  More  specifically,  if  we  let  Pe(f) 
represent  the  estimate  computed  by  either  the  in¬ 
direct  or  the  direct  approach,  we  have 

*vg[P#(f)j  =«  Q(f )*P(f ) 

where  the  average  may  be  taken  either  over  the 
ensemble  or  along  time  [ 3 J .  Q(f)  is  called  the 
spectral  window  corresponding  to  the  window  used 
in  the  computation. 

The  direct  and  indirect  approaches  to  esti¬ 
mating  the  power  spectrum  of  a  stochastic  process 
are  not  fundamentally  different.  Blackman  and 
Tukev  [31  have  derive!  the  ccv.cc t ic’  between  a 
data  window  Wft)  and  its  m valent"  lag  win¬ 
dow  b(~3.  this  relat i onsh i p  turns  out  to  be  the 
correlation  integral 

l.(~3  -  ~  Wfl)W(tJT)dt  . 

By  "equivalent"  we  mean  that  tbc.  appropriate  use 
of  the  respective  'Vquival  ;nt"  windows  yields 
power  spectral  estimates  which  are  equal  in  ex¬ 
pectation.  Note  that  if  the  data  window  is  an 
i-veii  function  of  time,  then  the  equivalent  lag 
window  is  proportional  to  the  convolution  of  the 
data  window  with  itself. 

At  this  point  it  should  be  clear  that  in  the 
indirect  method,  the  spectral  '.’indew  Q(f)  cor¬ 
responding  to  the  lag  window  Lf'O  is  simply  the 
Fourier  transform  of  For  the  direct  method 

it  then  follows  that  t?»e  spectra!  window  Q(f) 
corresponding  to  the  data  window  W(r  3  is  given 
bv  Q(C)  =  rp  I  I  f  t  )  )  *  where  Iff')  is  defined  to 
he  the  Fourier  transform  ol  W  ( t  3  .  J(f3  is 
called  the  frequency  window  correspondin'  to 
’rt(t).  We  again  refer  the  reader  to  Blackman  and 

Tukey  (31. 
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class  so  tlt.it  lh*'  rw'slifi 


3,  srnr.T  data  b ' bnov.N 
M  motivation 

lr  1  f  »’1i  t  of  the-  rat  ho  i  simple  i  el  at  ionshi ;> 
hot  Vi  on  a  data  window  and  its  equivalent  lav  win¬ 
dow,  one  should  not  ignore  the  i*>:t  ens  i  ve  ciffit 
which  has  boon  pul  into  the  construction  of  lag 
windows  hv  numerous  nut  hots  nor  the  properties  of 
Lhose  windows.  See*  for  exarpJe  the  cor;?i  Jr.  ■  ir.i 
given  by  Jenkins  f  9)  or  by  1'arzen  [12,13). 

For  instance,  recall  the  continuity  c ) a.*s  to 
which  these  functions  belong.  In  parti  culm  we* 
observe  that  one  of  the  lag  windows  which  has  been 
proposed  by  Parzcn  f 13]  is  proportional  to  a  di  - 
latcd,  fourth  order  B-splinc  bRsis  function.  One 
of  the  two  data  windows  suggest  oil  by  Welch  [18] 
lias  the. shape  of  a  second  order  r  Bp-spline  basis 
function;  he  recognized  that  his  proposed  window 
yields  the  Parr.cn  spectral  v:indow.  Henceforth 
we  will  refer  to  this  window  which  is  given  bv 


as  the  Parzcn  data  window,  even  though  it  is  also 
a  member  of  the  family  of  spline  windows  proposed 
in  tin’s  paper.  Today,  Welch’s  procedure  (de¬ 
scribed  in  the  introduction)  incorporal i ng  the 
Parzcn  window  is  undoubtedly  the  most  widely  used 
method  of  direct  power  spectral  estimation. 
Bingham  et  al  (2]  proposed  the  following  data 
window  which  is  continuous  through  the  first  de- 
rivat i ve: 


It  I  <  .AT 


W(t) 


.at  <  It  I  <  .yr 


N  0  ,  otherwise  . 

Wc  have  normalized  this  window  properly  although 
this  was  not  done  in  [2], 

From  another  point  of  view,  wc  note  that 
under  the  Fourier  transform  operator  the  effect 
of  an  occurrence  at  one  point  in  time  tends  to  be 
spread  over  all  frequencies,  Fur rhermoro .  time 
functions  with  corner!!  have  Fourier  transforms 
that  have  a  "ringing"  effect.  Thus  It  seems  de¬ 
sirable  to  have  a  data  window  of  a  high  continuity 


V’ 


class  so  that  ill*1  nodifio*  sample  function 
smoothly  to  /ir*i  at  it!:  nuls. 

Such  coos  i  cut  at  i  r.otivut.  ore  to  e/ar  :r-* 
the  p .  s  i  b  i  1  i  l  v  of  usi’v  dal  a  uii  ovs  b«v»cM  upon 
higlu-r  oriii  i  li-spliue  bPsit  i u:.<  i  i .  ^  .  A  prdin- 
Iniifj  survey  is  ev.rout  ay  i  up,.  1‘iist,  ve  ray  c  hoer» 
the  desired  continuity  class.  Stcondlv,  lh.  coi 
poi.diny  spevt  ral  windows  have  th<  shape  ^  Hi — J 

so  that  the  magnitude  of  the  s  i  Jo  h*  )>*■•>  rsjtidlv 
decreases  relative  to  the  rain  pi  nr:  h *:  the  con*  i  n* 
uity  (or  order  k)  increases.  It  addition,  all 
spectral  windows  are  of  ncn-ner.et  i ve  type  (Pitrzen 
f  1  5 1  ) .  l’hirdly,  deboot  [7]  lias  derived  a  genet  a  1 
algorithm  which  can  be  used  for  accurately  evalu¬ 
ating  any  order  li-spline  basir;  function  \.*itb  anv 
Interval, of  snppott  and  at.  any  nu-r.b  .length, 

3.2  TliH  rUNoAMKHiAl.  SH  I*::  IT’SH  IONS 

Schoenberg  (  1A  ]  introduced  the  spline  funr. 
tions  in  19  A  A  with  the  fol  lewi  ng  definition.  V.’ith 
k  a  positive  integer,  o  real  function  a/k(l)  de¬ 
fined  for  nil  til*'  is  called  a  spline  functii*. 
of  order  k  or  degree  k-1  if  it  has  th.'  folio.** 
ing  properties: 


(O 

(U) 

(Ui) 


®4(t)  Is  of  class  Ck  t(F). 

is  composed  of  polynomial  arcs  of 
degree  at  most  k-1  . 

The  polynomial  arcs  arc  joined  novheie  other 
than  possibly  at  integers  n  ff  k  is  even 
or  at  points  n4^  if  k  is  odd.  The 
spline  is  said  to  have  knots  at  these  points 
I f  we  define 


K-  • 
u 


y  >  0 


l-i 


y  <  0 

with  exception  of  the  special  case  0+  ~  £  and 

let  6k  represent  the  k*^  eider  central  ditte:- 
encc  operator  with  unity  step,  i.e.,  ft^ffy)  * 

f(y)  and  6kf(y)  -  6k~5  f(y+*)  -  6k"1f(y^) 
for  k  -  1,2,*..,  tlicn  the  function 

M,(t)  -  ~*)T  eti*1 

Is  easily  shown  to  be  a  simple  example  of  a  spline 

function  of  order  k.  M*(t)  is  called  the  funda- 


mental  B -spline  basis  function  of  order  k  .  We 
list  a  few  of  the  many  important  properties  of 
Mk(t)  and  refer  the  reader  to  Schoenberg  [  , 

Curry  and  Schoenberg  [6],  and  to  deBoor  [7j  for 

these  and  other  properties.  M  v  ( t )  is  positive* 

k  k 

in  the  open  interval  (-  )  and# ident ically 

zero  elsewhere;  it  has  knots  at  -  j+v  where  v 
is  an  integer  satisfying  0  <  v  <  k.  Mk(t)  is 
an  even  function,  and  its  integral  over  P  is 
unity.  Mk(t)  also  has  the  representation 


Mk(t)  * 


T^TyT^ 


-|+, 


r 


For  example,  the  fourth  order  function  is  given 
by 

r  i  /  .  o  r  .  .\  _  .  . 

Iti  <  i 


K*(t)  = 


< 


, l<  It  I  <2 


.0  , otherwise. 

We  shall  require  two  integral  propertied: 


M.<f)  -  j 


and 


CVM^Kt)  -  MpM(T)  . 

3.1  V\\K  PROPOSED  CLASS  OF  DATA  WINDOWS 
3.3.1  Ini roduc t i on . 

The  class  of  spline  data  windows  we  propose 
consists  of  the  fundamental  spline  functions  each 
dilated  so  that  Its  support  is  equal  to  the  length 
of  the  sample  function  and  normalized  so  that  the 
total  power  of  the  data  remains  invariant  when  the 
spline  data  windows  are  applied.  If  we  let  Ck 
iopreseat  Liu*  normalization  factor,  then  we  can 
use  the  results  of  .the  previous  sections  to  ob¬ 
tain  the  following  table  of  Ju.»  spline  data  w in¬ 
flows  and  at*  th.:  corresponding  frequency  windows, 
lag  windows,  and  spectral  windows. 

Since  the  total  power  >r  is  given  by 

nr/ 

-  *  x*  (t)dt  *  T  Var[x|, 

-T/? 

we  seek  constants  C*  depending  upon  the  order 
k  of  the  spline  data  window  such  that  also 


■  -V? 


tcAe-Mor'dt 


where  we  use  the  tilde  to  indicate  an  un -normal¬ 
ized  window.  The  expectation  ot  this  equation  is 

*  -  C;  i‘  ‘  ElU\(0]K[xJ(t)ldt 
‘-"/a 


-  C;  Var(x|  '  Q;(t)dt. 

V‘T4 

Changing  the  variable  of  integration  we  obtain 

/  ^ 

c;  =  k/f  M?(t)dt  . 

1  • 

■3 

An  application  of  Pars oval's  theorem  yields 


0?  -  rk/l 


du . 


Therefore,  we  can  calculate  the  appropriate  nor¬ 
malizing  factor  for  any  order  spline  data  window: 
Cj  «  1.0  and 


Ck 


2  V  L-iliOL-ii-itlll 

.  i : t2k-t-i) : 


for  k''-^ . 


Normalizing  for  invariance  of  total  power  is 
equivalent  to  the  requirement 

1  (0)  -  j’"5  Q(f)df  -  l  . 


Data  Window 
Wk(t) 

c*V*(t)  ■ r* ■'•(¥) 

Frequency 

W i ndow 

•f*(0 

CVT  f i i n( -i T /k)  \  * 

‘  k  '  \  ~n7k  ) 

Equ  i  valent 

Lag  Window 
U(r) 

T  (u1  V’*<T)] 

c?;  /k-\ 

v-(f) 

Spectral 

W  i  nd  ow 
fh(l) 

m-j-)  i 

i 

<1  /  .  i«K  -fT/k)  \  2,1 

TAB  If.  1. 
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We  point  out  th;it  is  simply  a  trunca¬ 

tion  function,  l.j  is  the  lag  window  suggested 
hy  Bartlett  fl],  il*'d  Q\  is  the  corresponding 
liar i  Lett  spectral  window.  W3,  La ,  and  Q.,  are 
respectively  the  Parzen  data,  lap.,  and  spectral 
windows , 

For  the  problem  of  balancing  i  lie  resolution 
and  stability  of  power  spectral  estimates,  ex¬ 
plicit  expressions  for  the  bias  and  variance  of 
smoothed  estimators  have  been  developed  in  the 
HtcraJture.  Jenkins  and  Watts  [10]  and  others 
give  these  expressions  for  the  indirect  computing 
method.  Welch  [17,18]  derives  the  expected  value 
ami  variance  of  the  discrete,  directly  computed 
estimates  for  a  stochastic  process  which  is 
it.Hiss i an,  in  addition  to  our  previous  assumptions 
about  the  process.  In  each  case  some  property  of 
i ho  window  being  used  is» expl ic i £ ly  involved, 
the  practitioner,  however,  needs  seme  device  to 
aid  in  the  so l ec t  i on  of  a  v i e.d.ow  that  will  viol d 
i  compromise  suitable  for  his  purposes.  Tut!  most 
cc.tnon  measure  of  windows  is  c  tiled  bandwidth  (by 
j  a.)  logy  with  filtering  p  rob  lens)  I 

Several  definitions  are  used  tor  bandwidth. 
Hackman  and  fukey  [1J  use  inst»  ad  the  term 
equivalent  width,  which  they  define  to  be  the 
ratio  of  the  square  of  the  integral  of  the  spec¬ 
tra  t  window  to  the  integral  of  its  square.  if 
we  define 

l  P  Wa(Hdr, 

V  -30 

tin  n  the  equivalent  width  is  given  by 

»,(<))  [[;"%(!)>'• 

leak  ins  anti  Watts  110}  define  Midwidth  to  be 
the  width  of  the  rectangular  (spectral)  window 
fur  which  the  variance  ot  the  iodirectlv  computed 
ni-»i  tril  estimate  is  tlie  same  .»■;  it  is  f  >r  the 
aiven  window.  I  hi  a  definition  of  bandwidth  \ield- 
I  I'M 

IS'  (Q)  \  ■ 

A;t  Ins  In-en  noted,  lor  a  properly  normalized 
window  WU)dt  l;  therefore  Iff  .  We 


mention  \\[  for  its  application  to  f lie  spline 
lag  windows  given  in  Inblc  1.  However,  its  de¬ 
fining  pi  ope  i  t  /  is  no  long*  i  tine  for  Llie  direct 
method  of  computation  which  we  advocate.  Kor 
the  class  of  spline  windows  ve  have 


Q?(f)df 


•o:kvr  -  i)?(wv*l> 

k"  i!i4k-l-l)* 

thus  Bj  (Q )  is  easily  determined  by  inverting 
lk  .  By  using  least -squares  technique*  for 
k  -  3 , , . . , 34 ,  we  obtain  the  approximation 
MQ*)  (l.b»(k)e*-t0-  O.V2)/T  . 

Parzen  (13J  defines  bandwidth  as  the  width  of  a 

rectangle  which  has  the  same  area  and  sane  n.axi- 

r/.u.n  height  as  the  glv-a  sp««-tr»l  wir.d*«.\  ll.*  rc- 

furo  '  *  ntf)df 

B.  *  0  >  - — 

umv  Q(f) 

f 

and  1  o r  the  spline  win. !,/.%*»  we  have 
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variance  is  invns*d  :  pi  oport  i  *mal  to  i  he  hand- 


.  dth.  It  can  also  be  arj;u  lhflt  the  bnndu  j  *1 1  it 

t  :  1 1 .  window  riisl  bo  at  le.v;t  as  si;  ill  as  1 1 .<> 
wi*'!!.  ui  the  1 1 '! t  l  owe  s t  iraport-nt  detail  ia  tfu 

*  t  i  ii-i. 

i  ab to  2  is  a  compilation  or  the  vaiious  cou¬ 
sin  its  assoc  iati-J  with  the  spline  windows  given 
in  1  with  order  up  t  a  k-21  . 


Order 

Gone  t  ui’ t 

i  Band 

;i  dtl.s  >  i 

K 

Ck 

!  !')'■  i  (Q„ ) 

■  }'v  «Q  . '•  n 

1 

1 .0^:0 

l.SOO't 

1 .  »ldd:'> 

2 

1 .  7  Ys  0 

1.8543 

1.3313 

1.2757 

3 

2.33‘.s 

!  2.301/ 

1 .  f  f  dd 

1 .  S  7  !  1 5 

4 

2.  its; 

2. 1.8  s; 

1.917 

1 . 8701 

5 

3.4093 

3.0203 

1  2  .  1  57  I 

7 . 0  <  9 

..  .  6 

3.9.V7 

3.3214 

1  2.  v.3e 

7  .  ?  •’  » 7 

7 

4 . 3  7  7  0 

3.4974 

2.4  IDS 

8 

4.8348 

3. 8  Sib 

;  2.7379 

7  .  Srf S3 

9 

5.2785 

4.0939 

2.9071 

.  7  4  3; 

ID 

S .  7 1 07 

4 . 37  09 

j 

3.0669 

1 

2.89  id 

11 

b.m? 

| 

4 . S3  6  5 

3.21^8 

t .  CM  M 

\? 

6.54 il 

4 .  j  4  2  i 

3.3t.3't 

1 .  1  /  d  t 

n 

6.97 3 8 

4.9396 

3.5030 

3.3011 

14 

7.3413 

S . 1292 

3.6367 

3.4264 

IS 

7.7293 

S.3121 

3.7657 

3 . 54 73 

lb  1 

8.1118 

5.4890 

3.8904 

3.61.41 

17 

8.4880 

5.6602 

4.0113 

3.7774 

18 

8. 8487 

5.8265 

4.1286 

3.8873 

19 

9.2243 

S. 9882 

4.2427 

3 . 9°4  3 

20 

9 .  SS  SI 

6. 14 S6 

4.3538 

4.0984 

?1 

9.9414 

6.2990  1 

4.4621 

— - 1 

4.2000 

TAB  IF.  2  . 

3.4  COMFIT  I  NO  WITH  STUNT.  DATA  WINDOWS 
Except  fot  the  choice  of  data  windows,  the 
method  we  propose  for  computing  power  spectra  is 
that  of  Welch  {181.  Departing  from  our  conven¬ 
tion,  we  assume  that  we  now  have  a  sample  func¬ 
tion  x(l)  ov ->r  the  interval  0*:t^S.  Using 
for  example  a  bandwidth  criterion,  choose  the 
order  k  of  the  spline  window  to  be  used  and 
then  fh<-  wind'"*  length  7.  Subdivide  y(t)  into 


N  segments  each  of  length  T  and  with  a  d i s - 
p!;.r<uK'iU  D  between  the  initial  points  of  cos- 
somtive  sep:;«  nt  s  .  Denote  tin*  Sty.  M.nts  by  >:  t  ( t  )  . 
Tims  (N  -  1)1;  i  V  -  S,  and 

x;  (t  )  --  x  -i  (  j  -  I )[)  -t  ~  j  for  j  1 ,  . .  .  ,N;  f :  .*  . 

Welch  hsis  the  choice  of  I)  onlv  as  a  paranoic* 
to  increase  the  si  ability  of  the  it  jolting  tsii- 
ne-.tid  power  spec  Li un,  in  the  caii  c> f  split.)  win- 
do.;*:  ;m  additional  < one idornt ion  should  be  made. 

!  1  I)  -  T ,  s  oiue  o  i  the  data  is  wo  i  i  ,h  t  ed  substan¬ 
tially  heavier  than  the  rennindet  .  As  can  be 
Si  i  n  from  Figure  1,  the  inequity  is  more  pro- 
no.  need  as  the  order  of  the  spline  window  in- 
crea.es  in  the  sense  that  the  spline  data  win¬ 
dows  become  prop.re*  slvoly  norc  peaked.  Since 
one  normlly  assr.e.-'  that  all  port  ions  of  his 
data  are  equally  representative  of  the  stochastic 
process,  intuitively  the  order  of  the  spline 
window  should  also  influence  I>.  In  practice  a 
pond  rule  of  thivb  for  relatively  short  sample 

functions  seetrs  to  b^  0  -  *T  for  V*  “  A  I)  ^  T 
1  * 
fui  k  -  10,  and  I)  ^  yT  for  V.  -  ,:0  with  the 

ovei  lap  for  other  order  window:’  la.^d  on  these. 

C 

Next  compute  the  smoothed  spectral  estimates 

F.(f)  -yl;  K)(f  )Mk  (t  }  “* .  j-1 . N. 

1  "  ~'fi 

Finally,  the  est itnat cd  power  spectrum  is  given  by 

r.(0  -  £  *  Mf>  . 

w  J-  1 

Regardless  of  the  order  of  spline  window 
used,  P,(f)  is  obviously  a  consistent  estima¬ 
tor  of  P(f)  (at  least)  whenever  the  stochastic 
process  is  Gaussian.  'lo  see  this,  the  express¬ 
ions  for  £(!*,(  1)1  and  Var[l*#(f)l  as  given  bv 
Welch  117,18]  need  only  be  compared  with  the 
convergence  properties  of  Liu*  spline  windows  as 
di ’.cussed  in  this  paper,  lb.  at  i ,  as  S  v*  ve 
have  the  following  type  convergence  for  a  sta¬ 
tionary  Gaussian  process:  if  T*-«a  such  that 

£-0,  then  E((P(f)-  F,(D)a}-0  . 

In  implementing  this  procedure  on  a  digital 
computer,  one  could  derive  a  closed  form  for  the 

R-':jlinr  basis  functions  like  that  pre- 
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viously  for  M^ft).  However,  for  a  general  pfo- 
grnm  and  ono  that  would  undoubtedly  bo  numeri¬ 
cally  preferable* ,  the  algorithm  of  deftoor  (7)  is 
recommended  for  evaluating  the  spline  functions. 

Jt  courso,  by  symmetry  only  halt  the  window  need 
be  directly  evaluated.  Brigham  (4]  thoroughly 
discusses  the  use  of  the  VV1'  for  compnt  ing  the 
discrete  Fourier  transform  which  approximates 
the  Fourier  integral  we  used  in  writing  Pj(f). 

4 .  EXPERT MENTA h  RKS U LT3 

The  effectiveness  of  the  spline  data  windows 
is  readily  apparent  from  experiments  with  machine 
generated  autoregressive  processes.  Power  spec¬ 
tral  estimates  computed  using  various  order  spline 
data  windows  are  compared  with  those  which  used 
the  Bingham -God f roy- Tukoy  and  Parzon  data  windows 
and  with  estimates  computed  by  the  indirect  method 
jf  Blackman  and  Tukoy  incorporating  the  Hanning 
lag  windi  j  [31.  The  estimates  were  nut  :;al  i/.ed  bv 
dividing  out  the  virl.rce  of  ?h«  toop  .iMc  pro¬ 
cess;  as  is  customary  when  one-sidid  power  spoetra 
a»*e  used,  the  magnitude  has  been  doubled. 

For  these  experiments  we  used  sample  func¬ 
tions  of  -ui to r ogres s ive  processes  which  were 
machine  generated  according  to 


covariance  fund  ion. 

(ii)  a-G-*-  :*m:  Hiugl.n.i -t.od f rey-Ti.ko/  d.ji  a 
Window  was  use. I  vi  th  segment. s  of  the 
•sa.-pJe  function  of  length  T  *:< 
samp  1 es . 

(iii>  Parzon  xx;  Parzen  data  window  was  used 
•with  segments  of  the  sample  function  of 
leu;;th  T  •  xx  samples. 

(iv)  tipline  xx-yv-zz;  Spline  data  window  of 
Older  k  —  xx  was  used  with  segments  of 
the  saw.plc  function  each  with  length 
T  -  yy  samples  and  displacement  D  -  zz 
samples  between  initial  points. 

Mte  graphs  illustrating  the^o  experiments  arc 
found  in  Appendix  1. 


at*.  t  1  ,  t  •*  l,...,ltb>4 


where  m  is  the  order  «»f  the  process,  (ol }  .ire 
the  regression  coefficients,  and  e,  is  a  white 
noise  with  Ef-\  1  U  and  Vnr(s,  )  -  1. 


i  in’  par titular  im  oc« 


.vaich  wi  re  used  in  the 


i  I  l  u  ;  l  ri*'  ions  had  the  following  coefficients: 

(i)  First  order,  -0.40 

(ii)  Second  order,  ..  I .00,  -0.30 

(iii)  Fourth  order,  v  0.7  3,  » „  0.30, 

,,  -  -0.33 ,  ...  *  -0.30  . 

Ihr  following  scheme  idr-uiifies  the  windows  in- 
vu!  e  <l  In  t  tm  computation  of  the.*  various  power 
spectral  »■  .t  i  mates. 

fi)  dinning  xx;  estimate  computed  by  the  in¬ 
direct  . tFv>*!  with  the  Hanning  lag  window 
[1]  and  using  x<  lags  of  the  simple  auto- 
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